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1. Introduction 

In 1983, Hawking proposed that, due to quantum gravitational effects, a scattering 
process might result in an initial pure state becoming a final mixed state [1]. This idea 
(which had been suggested earlier as a phenomenological possibility worthy of investiga- 
tion [2]) is clearly one of great intellectual power. It would represent a truly fundamental 
modification of the laws of physics. 

Nevertheless, Hawking's proposal quickly received criticism. It was noted [3,4] that 
conservation laws become decoupled from symmetry principles, so that (for example) ro- 
tation invariance would no longer imply conservation of angular momentum. Furthermore, 
actual calculations in models which admit systematic approximations showed no hint of 
the decay of purity [3] (see, however, ref. [5]). Also, since quantum mechanics is well tested 
in a variety of systems, corrections to the generalized hamiltonian would have to be quite 
small, with eigenvalues of order 10^ GeV or less [4]. 

By far the most severe attack was leveled by Banks, Peskin, and Susskind (BPS) [6], 
who concluded that conservation of energy and momentum could not be preserved without 
losing locality. Since then, however, the prevailing notion of the meaning of locality has 
weakened somewhat, due to the discovery of the wormhole phenomenon [7]. Wormholes 
violate locality, but the phenomenological results turn out to be remarkably benign [8- 
10]. It therefore seems appropriate to reexamine Hawking's proposal, to see whether or 
not the breakdown of locality noted by BPS is truly objectionable, and to see what other 
constraints, if any, may be placed on the decay of purity. 

Let us briefiy review some of the basic formalism for a theory in which pure states 
can evolve to mixed states; for more details, see refs. [2,4,6]. Such a theory must be 
formulated in terms of a density matrix p rather than a state \il)). The eigenvalues of p 
represent probabilities, and so must be real, positive or zero, and sum to one (assuming 
normalizable states). Thus p must be hermitian and satisfy Trp = 1. The state is said 
to be pure if p = \ip){ip\ for some Hilbert space vector \ip); in this case Trp = Trp = 1. 
If Trp < 1, then p cannot be written as \ip){ip\ for any IV'), and the state is said to be 
mixed. The usual evolution equation for p is p = —i[H^ p], and this leaves both Trp and 
Tr p constant in time. To allow loss of purity, we must modify this equation. Following 
refs. [2,4], we assume that the new evolution equation for p is still linear, and still first 
order in time derivatives: 

Pab = $ab^ Pdc- (1) 

The generalized hamiltonian ^ must be constrained to preserve the hermiticity, positivity, 
and trace of p. As shown by BPS, the most general equation preserving p^ = p and Tr p = 1 
is 

p = -i[H, p] - WpiQ^'Q^P + pg"g^ - 2Q^pQ") , (2) 



where the Q's are any hermitian operators other than the identity, and (/^./j is a hermitian 
matrix of couphng constants. The hamiltonian H is unambiguously defined as the operator 
appearing in the commutator term. BPS also showed that a suflficient condition for p to 
remain positive is that the eigenvalues of g^^p be positive or zero. 

Acceptable phenomenology leads us to demand that the familiar conservation laws be 
satisfied. Energy conservation, for example, requires that Tr f{H)p — 0, where f{H) is 
any smooth function of the hamiltonian H. BPS have shown that, H QaB has nonnegative 
eigenvalues and is, in addition, real and symmetric, then conservation of energy requires 
that H commute with each of the Q's. (If g^,/] does not satisfy all these conditions, there 
may be other possibilities; this will be discussed later.) In a quantum field theory, there 
are very few operators available which commute with the hamiltonian, and all of them are 
global: that is, integrals over all space of a local density.* We always have at our disposal 
the hamiltonian H and the total momentum operator P to press into service as Q's. We 
may also have global charges corresponding to conserved quantities like baryon number. 
Is eq. (2), with the Q's chosen from this list, a viable possibility? 

BPS concluded that the answer is no, due to a breakdown of locality. In particular, 
they state that cluster decomposition will no longer hold. Given the density matrix p at 
some fixed time, and localized operators A{x) and B{y), cluster decomposition states that 

Tr[pA{x)B{y)\ ~ Tx[pA{x)\ Tr[pS(y)] (3) 

when \x — y\ becomes large. The precise meaning of "large" depends on the theory, and on 
p. In ordinary field theory, eq. (3) is usually quoted as a property of the vacuum density 
matrix pg = |0)(0|. If the lightest particle in the theory has mass to, and if we take p = pQ, 
then "large" means much bigger than mT . But if, for example, p represents a pure state 
consisting of one or more particles whose wave functions are localized near the origin but 
coherent over a region of linear size a ^ mr\ then "large" means much bigger than a. 
The question is, if eq. (3) is valid for some initial density matrix pinit, and for \x — y\ 
greater than some length L, will this still be true at later times? 

If we assume that energy is conserved, then pure energy eigenstates like po do not 
evolve in time, and so there is no modification of the cluster property for p-^^i^ = Po- If 
we consider instead a pj^it corresponding to a localized excited state, we must ask how 
quickly its wave packet spreads out as time evolves. Even in ordinary field theory, wave 
packets spread out, and the minimum value of \x — y\ for which eq. (3) holds will grow with 
time. It turns out that eq. (2) does cause wave packets to spread out a bit faster than they 
otherwise would, but the effect is not dramatic; certainly it does not qualify as a violation 
of cluster decomposition. However, if we put some particles at both x and y, then extra 
correlations do develop with time. For free nonrelativistic particles, these remain small at 



In a gauge theory, we have the local generators of the gauge symmetry; these annihilate all the physical 
states, however, and so would vanish in eq. (2). 



all times, and actually disappear at late times. This is discussed in detail in sect. 2. We 
conclude that the nonlocality noted by BPS, while present, is in fact difficult to detect, 
and that eq. (2) cannot be dismissed because of it. 

One may still question whether or not eq. (2), with global Q's, has any reasonable 
chance to arise as the low energy limit of a more fundamental theory. I know of no such 
theory, but there is a strong passing resemblance between eq. (2) and the effective theories 
that arise due to wormholes [7-10] (which were originally thought to cause purity to decay). 
In wormhole theory, the usual field theory action Sq = J d x C is modified to [10] 

^ = ^0 + QapQ^'Q^ + ■■■ (4) 

where Q^ = J d x q^{x)j q^{x) is some local operator, and the ellipses stand for higher 
powers of Q's. It seems not totally unreasonable to suppose that there could exist a variant 
of wormhole theory which ultimately results in eq. (2) rather than eq. (4). Of course, there 
are some important differences between eq. (2) and eq. (4): in eq. (4), the Q's are four 
dimensional integrals, and they modify the action itself. They therefore do not lead to 
violations of quantum mechanics, but do look like they would produce severe violations of 
locality. As is now well known, this is not the case. Similarly, the apparent violations of 
locality in eq. (2) are not as drastic as they appear to be at first glance. 

If eq. (2) does not lead to serious violations of locality, and cannot be immediately re- 
jected on grounds of implausibility, it becomes necessary to reexamine its viability. That is 
the purpose of this paper. Phenomenological constraints, including the anomalous spread- 
ing of wave packets, are examined in sect. 2. Formal restrictions are discussed in sect. 3, 
where it is found that eq. (2) is incompatible with the weakest possible form of Lorentz 
covariance (if energy is conserved). This is the main conclusion of this paper. However, a 
possible loophole in the argument is presented in sect. 4, and conclusions in sect. 5. 

2. Phenomenological constraints 

As noted in the introduction, the question of the locality of eq. (2) requires us to study 
the spreading of wave packets with time. Rather than immediately attacking this problem 
within quantum field theory, we will warm up with some examples from one-dimensional, 
nonrelativistic quantum mechanics that contain the essential physics. Let us first consider 
a harmonic oscillator of unit mass and frequency (so that the hamiltonian is H = a'a + ^) 
which is in a coherent state. At time t = 0, the state is specified by an arbitrary complex 
parameter l3 = |/?|e*'^: 

|/3,0)=e-|^l'/V«^|0). (5) 

It is well known that at later times the state is still coherent, but with a different value of 

(3: 

|/3,t) = e-^^/2|/3e-^0). (6) 



Furthermore, the probabihty to find the osciUator at position x at time t is given by 

{x\p{(3,t)\x) = 1 e'[^-^^l/5|cos(t-0)]' (7) 

where p{(3,t) = |/3,t)(/3,t| is the density matrix at time t. Thus, the oscillator's wave 
packet just moves back and forth without changing shape, its center following a classical 
trajectory with maximum amplitude y^ |/9|- 

Now consider modifying the quantum equation of motion to the general form of eq. (2). 
In order to conserve energy, the Q's must commute with H. For simplicity, we will take a 
single Q and set it equal to H, and let the corresponding value of g be real and positive. 
To compute the density matrix at time t, we first expand the initial state \(3, 0) into energy 
eigenstates: 

cxj on 

|/?,0) = e-|/^lV2^^|,), (8) 

where \n) = (n!)^ ' (a')'^|0) is a normalized energy eigenstate with energy n + ^. The 
time evolution of a density matrix element \n){m\ in this basis is simple, and results in 



oo 



n o*m 



(3^(3 



The probability of finding the oscillator at position x at time t is again (a;|p(/3, t)|a;); let 
us study the behavior of this function as t ^ oo. Looking at p{(3,t), we see that its 
off-diagonal components decay away at late times, and it approaches 

p(/?,oo)=e-|^l J2-^\n){n\. (10) 

n=0 ""■ 

Furthermore, this final state is independent of the specific choice of the Q's as long as 
they all commute with H. To compute {x\p{(3, oo)|a:) explicitly, note that the off-diagonal 
components of p{(3, t) can also be removed by averaging over the phase of (3: 



1 /'^'^ 
p(/3,oo) = — / #p(/3,t). (ii; 

^TT Jo 



Since the result is independent of time, we may as well use p(/3, 0) = |/9,0)(/3,0| on the 
right hand side. Then, taking the expectation value in a position eigenstate and using 
eq. (7) at t = 0, we find 

{x\piP.oo)\x) = -1^ /"'"c^^^e-t^-^l^l-^^]' . (12) 

27r'^/^ Jo 
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We thus see that the osciUator forgets its starting point, but conservation of energy forces 
it to remember its maximum amphtude \/2 \(3\. If we make a change of integration variable 
in eq. (12), we get 



{x\p{P,oo)\x) 



+V2I/3 

V2|/3| 



dy 






^^m'-y' 



(13) 



where appropriate factors of U have been restored. Here vr" (2|/3| — y )^ ' represents 
the probabihty for a classical oscillator with a specified maximum amplitude v2|/3| to 
be found at position y at a randomly selected time. Eq. (13) tells us that this classical 
probability distribution is then smeared out over the width of the quantum wave packet. 
In the h ^ limit, the wave packet becomes a delta function, and we recover the classical 
probability. 

This is a rather innocuous result. In order to discover this sort of violation of quantum 
mechanics experimentally, one would have to let the oscillator run for a long time (since 
the new parameter g is presumably small) while keeping precise track of the time. On 
the other hand, systems such as the earth revolving around the sun can be treated in an 
essentially similar way, have been around for billions of years, and have not been seen to 
exhibit any mysterious timing discrepancies. But it is hard to translate this into a limit 
on g for elementary particles in any convincing way. 

As another example, let us consider a free, nonrelativistic particle with unit mass in one 
dimension (so that the hamiltonian is H = ^p )■ (For a related discussion, see ref. [11].) 
Take the initial state to be a gaussian wave packet of width a and and mean velocity pq: 



iIj{x,0) 



1 



l/4„l/2 



exp 



n 



X 



2 n 



IPQX 



2o? 



(14) 



Then, at later times, we find (according to ordinary quantum mechanics) that the proba- 
bility to find the particle at position x is 



\i^{x,t)[ 



1 



^1/2(^2+^2/^2)1/2 



exp 



0? + t^/a^ 



(15) 



This is easy to understand. The initial uncertainty in position Aa; ~ a implies an uncer- 
tainty in momentum Ap ~ 1/a. Then, at later times, the momentum uncertainty implies 
an additional uncertainty in position (Ap)t ~ t/a. These are then compounded quadrati- 
cally to give the time dependent position uncertainty Ax(t) ~ (a + t /a ) ' implied by 
eq. (15). Meanwhile, the center of the wave packet moves with constant velocity pQ, the 
only possibility consistent with conservation of momentum. 

Now we again consider modifying the quantum equation of motion. Since the hamil- 
tonian is so simple in this case, we can take the Q's to be arbitrary functions of p and 
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still conserve energy and momentum. Let us work out one example in detail: a single Q, 
this time (for later mathematical convenience) equal to the momentum p. Expanding in 
momentum eigenstates, we see that the probability to find the particle at position x at 
time t is 



{x\p{t)\x) = — 



dp dk ij{p, 0)ij*{k, 0) e^b-fc)^-*(p'-fc')V2-gb-fc) V2 



(16) 



where i(;{p, 0) is the Fourier transform of '(/'(x, 0): 

al/2 



V'(P,0) 



n 



1/4 



exp 



^a'ip-PoY 



(17) 



Because the energies are continuous in this case, we cannot simply discard the off-diagonal 
terms as we did for the harmonic oscillator. Instead, we note that the choice of Q = p 
has rendered {x\p{t)\x) as a gaussian integral which can be evaluated exactly at arbitrary 
times: 



{x\p{t)\x) 



1 



exp 



(x-poty 



a^ + t^/o? + 25ft 



(18) 



7rV2(a2 + t2/a2+2(7t)V2 

We see that there is a new contribution to the uncertainty in position, Axnew ~ (2(7t) ' , 
which is compounded with the ordinary quantum uncertainties. The t ' dependence of 
Axnew is characteristic of a one-dimensional random walk. Note that, since momentum is 
still conserved, the center of the wave packet still moves at constant velocity pQ. 

This change in the nature of how wave packets spread out is again an innocuous one. 
In fact, for g <til, the effect of Axnew is negligible, since then gt <ti a -\-t /a at all times. 

Now let us consider two free particles, with hamiltonian H = ^{pi + P2) ^^^ initial 
wave function 



'(/'(xi, X2, 0) ~ exp 



ipiQXi + ip20X2 



(xi-xio)^ (x2-X2o)^ 



2a' 



2a' 



(19) 



We can analyze this system in the same way as before; taking Q to be equal to the total 
momentum pi + p2 of the two particles, we eventually find 



(xi,X2|/o(t)|xi,X2) ~exp 



^2/^2 



(g^ + tyg^ + 2gt) {yf + y|) - 4gt yiy2 
(a2 + t^/a'^){a^ + t2/a2 + 45ft) 



(20) 



where y^ = Xi — [xiq +Piot). The yiy2 term indicates that the presence of a second particle 
infiuences the first. Even though they are spatially uncorrelated initially, correlations 
develop, independent of how far apart they are. This represents a violation of cluster 
decomposition. However, for (7^1, the new correlations are small at all times, and 
actually decrease with time for t > a . Furthermore, if we trace over the states of the 



second particle (integrate eq. (20) over X2), we reproduce eq. (18). Thus the second particle 
has no effect if we know nothing of its state. 

A quantum field theory will exhibit similar behavior. Indeed, for free field theory this 
analysis can be taken over completely, except that we should use a relativistic dispersion 
relation, and choose Q's which respect Lorentz covariance. The change in the dispersion 
relation prevents exact evaluation of the integrals, and even a steepest descent analysis is 
quite involved, but the behavior should be qualitatively similar. This means that loss of 
purity does not force us to rethink the LSZ reduction paradigm for scattering theory, since 
wave packets will remain sufficiently compact, and extra correlations will always be small. 
A more serious problem is the requirement of Lorentz covariance, which, it turns out, is 
impossible to satisfy. This is the subject of the next section. 

3. Formal Constraints 

Since we have a theory in which correlations can appear between spacelike separated 
systems, finding a Lorentz covariant version may seem to be a hopeless task. This is 
not necessarily the case: theories with tachyons, whatever their other faults, are Lorentz 
covariant, and have superluminal correlations. Thus a reasonable question to ask is, what 
is the weakest possible form of Lorentz covariance which can be demanded of eq. (2) ? 

Let us begin with a concrete example of eq. (2) which is Lorentz covariant. Consider a 
scalar field ip with lagrangian C = ^d^^ipd^ip — V{ip) + Jip, where J is a gaussian random 
variable with {J{x)J{y)) = g 6 {x — y). Here angle brackets denote an average over the 
probability distribution for J. This theory is manifestly Lorentz invariant, and, as shown 
by BPS, leads to the evolution equation 

p = -i[H,p]-lg fd\[ipix),[ipix),p]] , (21) 

which is a specific form of eq. (2). This theory is not viable, since BPS demonstrated that 
it produces enormous violations of energy and momentum conservation. It does, however, 
serve as an example of the form eq. (2) could take in a Lorentz covariant framework. We 
can rewrite eq. (21) as 

aOp = -z[pO,p]-l<7y"dEO[v'(^),[v^(f),p]] . (22) 

Here x specifies a point on the spacelike hypersurface E whose unit normal defines the 
time direction, and dTi^ = -^el^^^ dx^dx jdx j^. Clearly, each term in eq. (22) carries a vector 
index in the time direction, and no other explicit indices. I wish to argue that this must 
be true of any version of eq. (2) which can arise from an underlying, Lorentz invariant 
framework. It is by no means clear that this condition is sufficient to ensure the existence 



of a full set of Lorentz transformations with the correct algebra, but it would appear to 
be necessary. This condition on eq. (2) is what I mean by the "weakest possible" form of 
Lorentz covariance. 

We do not need to know much more about Lorentz transformations to reach a strong 
conclusion. This is fortunate: since eq. (2) tells us that infinitesimal time translations are 
not generated by the hamiltonian H, it is clearly unsafe to assume that space translations, 
rotations, and boosts will still be generated by the usual operators.* If Lorentz invariance 
is to have any meaning at all, however, we must demand that the Lorentz generators 
(whatever form they may take) satisfy certain basic properties, including the following: 

(1) The total energy and momentum operators H and P should transform as a four- 
vector under infinitesimal Lorentz transformations (so that, for example, the generator of 
boosts in the ^-direction, acting on H, results in Pz). 

(2) The infinitesimal time and space translations dt and dx must also transform as a 
four-vector. 

(3) The usual product rules should hold, so that, for example, M = H — P is still 
invariant under Lorentz transformations. 

This information alone is enough to allow us to conclude that p and [H, p] both trans- 
form in the same way. I have argued that Lorentz covariance will be lost if the final 
term in eq. (2) does not also have this property. Therefore the product of two Q's must 
transform like H. The simplest possibility is to let one Q be a Lorentz scalar, perhaps 
M = H — -P , and let the other Q he H itself. The problem with this is that it implies 
a matrix g^p which is purely off-diagonal, and which therefore has a negative eigenvalue. 
One can see the problem explicitly in free field theory. For definiteness, take Q\ = H and 
Q2 = M , with gii = (722 = and gi2 = (721 = 9- AH one-particle states have the same 
value of M , and so the extra term in eq. (2) vanishes for all one particles states, pure or 
mixed. If, however, we consider a two-particle state \p1p2) = fl' (pi)«' (P2)|0), and evolve 
the initial density matrix p(0) = \piP2){P3P'i\j we find at later times 

p{t) = exp[-z(£;i2 - EM)t - g{Ei2 - Em){mI^ - Ml^)t\ \piP2){p?,Pa\ (23) 

where Eij = E^ + Ej, E^ = {pf + m^)^/^, Mf- = {Ei + Ej)'^ - {pi + pj)"^, and m is the 
mass of one particle. The g dependent term is not always negative; indeed its sign is frame 
dependent. (The same is true if we replace M by some other Lorentz scalar, such as a 
charge.) An off-diagonal matrix element of p which is monotonically increasing guarantees 
that p will eventually have a negative eigenvalue. Therefore this Lorentz covariant scheme 
must be rejected. 



In section 2, we implicitly assumed that the total momentum operator P generates space translations; 
relaxing this may provide another way out of unpleasant conclusions regarding nonlocal effects. 
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The only other possibihty arises in a super symmetric theory: then we have the super- 
charges at our disposal. They commute with the hamiltonian, and their product has the 
correct transformation properties. The appropriate version of eq. (2) is 

p^_,[H,p]-\g\Hp + pH-2[a''f''{Q{pQ^ + Q^pQ{)\ . (24) 



Trp = follows from the supersymmetry algebra {Qq,-:Q\} = (o"^) P^- (We may use 
conventions in which both (c )^^ and (a ) are equal to 5^0^.) Eq. (24) also conserves 
energy and momentum, since Q^ and Q^ commute with H and P. 

The problem with eq. (24) is that it does not conserve angular momentum. Although it 
is straightforward to show (using the transformation properties of Q^ and Q]^ that while 
Tr J p remains constant, Tr J p does not. It is easy to see this physically: Q ~ 6' a, where 
W creates a fermion and a destroys a boson. An initial state of one boson at rest will be 
converted to a state of one fermion at rest by eq. (24). This fermion is equally likely to 
have spin up or spin down, and so Tr Jp remains zero, but clearly Tr J p has changed. 
Thus this scheme must be rejected as well, and we are forced to conclude that eq. (2) is 
incompatible with the weakest possible form of Lorentz covariance. 

This objection would not arise in two dimensions, where there is no angular momentum. 
Perhaps eq. (24) has some role to play on the superstring worldsheet. 

4. A possible loophole 

The conclusion just reached — that the evolution of pure to mixed states is incompatible 
with the weakest possible form of Lorentz covariance — was based on the assumption that 
the Q's in eq. (2) must commute with H. BPS showed that this is the only possibility if 
g^Q is a real symmetric matrix with nonnegative eigenvalues. We must still consider the 
possibility that g^^jj does not obey these conditions (since positivity alone is a sufficient 
condition to avoid development of negative probabilities, and even that may not be neces- 
sary) . Is it possible to conserve energy, momentum, and angular momentum, even if the 
corresponding operators do not commute with the Q's? Surprisingly, the answer is yes. 
The example I have found is not a viable theory because it allows p to develop negative 
eigenvalues; however, I have been unable to prove that a model without this fatal flaw 
does not exist. This leaves open a small chance to maintain the weak form of Lorentz 
covariance along with the familiar conservation laws. 

The simplest theory which conserves energy even though [H, Q] 7^ is based on the 
hamiltonian H = h%^ where 6 is a fermion operator obeying the usual anticommutation 
relations {h,h} = and {b\b} = 1. We take Qi = ¥ + b and Q2 = i{b^ — b), with 
gil = —(722 = 9 and gi2 = (721 = 0. Obviously, (/q,^ has a dangerous negative eigenvalue. 
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Eq. (2) becomes 

p = -i[H, p] - 2g{h^ph^ + hph) . (25) 

We can now see why energy is conserved; if we compute Tr H^p, all the extra terms vanish 
because we always get a factor of either 6 = or b' = 0. However, it is not hard to 
show that, in general, an initially pure p immediately evolves into one with a negative 
eigenvalue.* 

The model may be made less trivial by extending it to a set of independent operators 

b^, with H = Z^blh■. 

p = -i[H,p]-2gj:{blpbl+b^pb;) . (26) 






It is easy to check that we still have Tr H^p = 0. This version can be promoted to a field 
theory of free (Majorana) fermions: 

p=-i[H,p]-2gj:jd^p[bHp,s)pb\p,s)+h.c.] , (27) 

where 6' (p, s) creates a fermion with momentum p and spin s. The extra term has the 
correct Lorentz transformation properties, but is nothing very simple (or local) in position 
space. And, as before, it leads to negative probabilities. 

This existence of this example is quite annoying, since it prevents a general proof of 
the necessity of [H, Q] = for energy conservation, yet is both untenable (because of the 
negative probabilities) and extremely contrived. It is hard to imagine a physical process 
(such as wormholes) leading to anything as complicated in position space as eq. (27). 



5. Conclusions 

A reexamination of the question of whether or not pure states can evolve into mixed 
states shows that phenomenological problems associated with loss of locality are less se- 
vere than previously believed. However, it seems impossible to maintain both Lorentz 
covariance of the modified evolution equation and the familiar conservation laws. A loop- 
hole in the general argument leading to this conclusion was noted, but the only example 
I could find which sneaks through it is rather special, and in any case results in negative 
probabilities. The fascinating possibility that purity may not be eternal is still out of reach. 

I would like to thank two anonymous referees for lively correspondence and helpful 
suggestions which substantially improved this paper. This work was supported in part by 
NSF Grant No. PHY-86-14185. 



After this paper was circulated as a preprint, Jun Liu found a modified version of cq. (25) which preserves 
the positivity of p [12]. 
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